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• individual islands and a statistical analysis of island shapes were used to reconstruct the evolution 

' of the island shape as a function of size. As they grow, islands undergo a transition from completely 

, unfacetted rough mounds (prepyramids) to partially {105} facetted islands and then they gradually 

evolve to {105} facetted pyramids. The results are in good agreement with the predictions of a 
recently proposed theoretical model. 



The morphology of the first three-dimensional islands appearing during strained growth of SiGe 
alloys on Si(001) was investigated by scanning tunneling microscopy. High resolution images of 
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PACS numbers: 68.35.Bs: 68.37.Ef; 68.55.-a 



|- ' The first stages of spontaneous island formation and evolution during growth of SiGe on Si(001) have recently 
\Q , attracted much interest i2i2i2i£i£iL£*2*i2i At typical growth temperatures, the Sii_ 2: Ge 2: /Si(001) system displays a 
sequence of morphologies qualitatively independent of x, at least for x > 0.2—i In the case of low misfit strain 
(x < 0.6), individual islands have been observed to form without nucleation 2,3 from a quasi-periodic array of ripples 
as a result of a strain driven instability>i2iiiiiii^ Islands appear as unfacetted mounds ( "prepyramids"—) on a stressed 
wetting layer (WL) and then they evolve into facetted pyramids bou nded by {105} facets as their size increases. 

Recently we have shown that the observations reported in Refs. I2I K1 ll l5l can be interpreted by assuming that the 
surface energy anisotropy of a SiGe layer on Si(001) has the following behavior: (i) orientations in a neighborhood 
. of (001) are permitted so that (001) is not a true facet, (ii) {105} is a facet orientation, and (iii) there is a range of 
unstable orientations separating (001) and {105}— With these hypotheses, the equilibrium shape of an island (in two 
dimensions) was calculated as a function of size. To simplify the calculations, {105} was not treated as a true facet 
O ■ in our model, although it is one in the real system. The important point, which is captured in the simplified model, 
is that there are a range of unstable orientations separating (001) and {105}. 

Depending on the island volume V there are shapes corresponding to prepyramids and/or truncated pyramids with 
a rounded top ( "T-pyramids" ) . The model predicts that there are only prepyramid islands for V < V\ , prepyramids 
^ . and/or T-pyramids for V\ < V < V3, and only T-pyramids for V > V3. 

Most morphological studies of the first islands appearing during SiGe growth on Si(001) have employed techniques 
with low spatial resolution, such as atomic force microscopy (AFM)£i&or low energy electron microscopy (LEEM)—^ 
The two exceptions studied pure Ge^ where prepyramid islands are too small to resolve much detail, or relatively 
dilute (25%) Gei& at relatively low temperature (for this composition), where the behavior may be rather different. 9 
We present here a scanning tunneling microscopy (STM) study of the pyramid precursors observed during growth 
of Sii_ x Ge x alloys on Si(001) in a composition range between x = 0.3 and 0.5. A detailed statistical analysis of the 
shape as a function of island size is used to verify the predictions of the above model. While a direct comparison of the 
data from three-dimensional (3D) islands to the predictions for two-dimensional islands is not possible, a qualitative 
comparison shows the striking similarities. 

The samples for this study were grown by ultrahigh vacuum (UHV) magnetron sputter epitaxy and investigated 
by room temperature (RT) UHV-STM. Heavily doped Si(001) substrates with a nominal miscut < 0.05° were flash 
cleaned by alternating current heating in order to remove the native oxide. A Si buffer 100 nm thick was grown to 
achieve a clean surface. SiGe layers with different composition were obtained by tuning the relative power applied to 
two sputter guns operating simultaneously. SiGe was grown at a substrate temperature T s = 600° C and a rate of 
about 0.08 nm/s. In some cases a short period of annealing at this temperature preceded the cooling of the sample 
to RT (at a rate of about 2°C/s) for STM characterization. 

Figure shows an STM image of a Sio.7Geo.3 layer 3.8 nm thick annealed for 20 s at T s . In spite of the very 



2 



different growth technique employed here, our results are compatible with those obtained by Sutter and Lagally 3 and 
Tromp et alSk who used chemical vapor deposition (CVD) and studied the film morphology by means of LEEM. In 
fact, Fourier transforms of STM images of this sample show that islands start to develop from a ripple structure 
composed of cells with an average spacing of 125 nm (see inset of Fig. [TJ. This finding is in good agreement with the 
value of 150 nm found in Ref. [3 for x — 0.3, within the experimental uncertainties of alloy composition. 

At higher resolution [Fig.Efa)] the mounds appear rounded and not regularly stepped, suggesting that, in contrast 
to clean Si(001), the (001) plane is not a facet for SiGe on Si(001). The contact angle between these prepyramids and 
the WL is zero, so it is difficult to identify the edge of the island base. 

When the Ge fraction x is increased to 0.5, the ripple structure is less clear, but still observable. As in the case 
of lower misfit, islands appear as unfacetted mounds [Fig. Hfb)]. As the composition is varied, prepyramids have 
qualitatively the same shape, but they occur at sizes decreasing with increasing x. Small prepyramids were in fact 
observed also for nominally pure Ge grown on Si(001) by Vailionis et al£ 

According to the model proposed in Refs.Q] and0, the cross section of small prepyramids is well fitted by a cosine 
function. The two islands shown in Fig.Elare close to the transition to the T-pyramid shape since small {105} facets 
are observable on their surface at higher resolution^ For this reason the line scans deviate somewhat from the cosine 
fits. We observe that the WL and the surface of the islands show the sa me patched (M x N) surface reconstruction 
with AI and N increasing with decreasing discussed in for the case of lower misfit alloys. 

In order to obtain insight into the shape evolution, we examined samples featuring islands at different stages of 
evolution (prepyramids, T-pyramids and mature pyramids 20 ) and measured their shape as a function of size. Samples 
of this kind were obtained by suitably choosing the amount 8 of deposited SiGe, having fixed the growth temperature 
to T s = 600°C and composition to x=0.5. 

At 8 = 1.5 nm the only 3D islands present on the surface are small prepyramids. At 8 = 1.8 nm a few facetted 
islands are observable together with partially facetted and unfacetted ones. However, a similar island distribution 
can be obtained after annealing a sample with 8 = 1.5 nm for 20 s at T s [see Fig. E^a) and (c)]. This observation 
is compatible with that of Jesson et al&Q and indicates that the critical thickness for the onset of 3D growth is 
kinetically determined, and that in general the whole film is not in equilibrium. The size distributions of the islands 
present in the above samples and in a sample with 8 = 1.9 nm were found to be compatible. This whole set of samples 
will be referred to as Sid, to denote that their surface is characterized by a relatively low density of islands (about 
10 9 cm -2 ). At later stages of growth the island density increases. Figure E2b) shows the surface of a sample (Shd) 
obtained by growing 2.8 nm of Sio.sGeo.s at 600°C and annealing it at this temperature for 20 s. Now, there are about 
10 10 islands per cm 2 . Besides mature pyramids (including islands close to the dome transition), small T-pyramids are 
observable also in this sample, as seen in Fig. [3d). 

Figure^Ia) shows two islands at different stages of evolution and gives an indication on how a prepyramid transforms 
into a T-pyramid. The upper island displays a small (105) facet appearing on regions of its surface of the with the 
greatest slope, as revealed by the cross-sectional line scan shown in Fig. Efb). Both the region where the facet meets 
the rough top of the island and that where it meets the island base appear curved, as predicted by the model discussed 
in Ref. The facet appears to sit on a shallow shoulder^ which shrinks for larger islands. This feature is in fact 
practically unobservable for the T-pyramid at the right lower corner of Fig.^Ja), consistent with the model. 

In order to describe the shape evolution as a function of island size, for each island we determined several parameters 
by means of a dedicated software program. As noted above, the bases of prepyramids are not sharply delineated. 
Roughness of the surrounding WL and "hut pits"^&24 [visible as small dark spots in Fig. EJc), and in more detail 
in Fig. 0] make it difficult to identify the island edge. For this reason, the island base perimeters were defined by 
means of a contour plot at a reference height ho- ho was chosen 0.5 nm above the average height at which {105} 
facets belonging to pyramids meet the WL. In this way the shallow shoulder (see Fig. 0} at the base of T-pyramid 
was not included in the analysis. Island heights were measured with respect to ho and volumes were calculated by 
integration. Two further parameters were measured for each island: (i) the area ^4(io5) of the island surface occupied 
by {105} facets (as described in Ref. l25j) . (ii) the area ^4(ooi) 01 the regions of the island surface with inclination less 
than 6 degrees. For prepyramids and T-pyramids, the latter parameter gives a measure of the total area and of the 
area of the rounded top, respectively. Islands that were too asymmetric, coalescing or poorly defined were excluded 
from the analysis in order to reduce the measurement uncertainty. 

We note that pits in the WL locally modify the strain configuration, possibly leading to the occurrence of asymmetric 
islands. At present we are not able to exclude that these pits are correlated to impurities 2 ^ introduced at the SiGe/Si 
interface during the necessary growth interruption before initiation of SiGe growth. If this is the case, it may be 
possible to produce samples with a flatter WL, allowing the uncertainties on the size measurement of prepyramids 
to be reduced and the model assumptions to be better approached. Finally the cooling of the sample may affect the 
shape of the islands and change it from its ideal equilibrium configuration^ 

Figure [S] illustrates the transition of prepyramids to T-pyramids and their gradual evolution to mature pyramids^ 
with increasing size. Prepyramids are characterized by a large A( oi) and a negligible A( 10 5), while the surface of 
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T-pyramids contains both, {105} facets and areas with orientation close to (001). 

Focusing initially on the samples S^, Figs.|3fa) and (b) show that the smallest 3D islands are prepyramids and that 
these islands "survive" in a limited range of volumes. The transition occurs at a volume of about 1000 nm 3 , where 
{105} facets appear [Fig.[3fb)]. Figure|5fa) suggests that a discontinuous drop in ^4( oi) occurs at the transition. As 
the island volume increases beyond the critical size, we observe that: (i) the rounded top of T-pyramids smoothly 
shrinks, so that the larger the T-pyramid, the sharper is its apex [Fig.OJa)] and (ii) the area of {105} facets increases 
monotonically as nearly a power law [Fig.JSfb)]. 

Qualitatively, the same considerations apply to the plots of A( oi) and A( 10 5) vs height [Figs.|S{d) and (e)]. However, 
while the range of volumes of prepyramids overlaps that of T-pyramids [Fig. [of a)], there is a clear separation between 
the range of heights of prepyramids and of T-pyramids [Fig. Eld)] . 

Thus from the data the following scenario is suggested. At small volume an island is a prepyramid. At some critical 
volume the prepyramid undergoes an abrupt transition to a T-pyramid, accompanied by a discontinuous jump in the 
area of {105} regions and a discontinous jump in the height of the island. Upon further growth, the island reduces 
the amount of (001) area while increasing the amount of {105} regions. 

Turning now to the data from sample Shd, Figs. Ufa) and (b) show that these islands have size distributions 
substantially different from those of samples Sid- The smallest islands observed in sample Shd [such as that pointed 
at by the arrow in Fig. Ofd)] are T-pyramids with a relatively large flat top and with a volume appreciably smaller 
than 1000 nm 3 . No unfacetted islands were seen in sample Shd, down to the smallest size that we can reliably identify. 
The model presented in Ref. |9| offers an interesting interpretation of this observation, as discussed below. 

The theoretical island shapes from Ref. were analyzed to determine shape parameters analogous to those measured 
in Fig. [S] We determined the "apparent" width of prepyramid and T-pyramid theory shapes from the points of 
maximum concave curvature [see schematic in Fig.|Sfa)]. For T-pyramids, the widths of the (001) and {105} regions 
were determined from the point of maximum convex curvature. 

Figures HJb) and (c) show the width of the (001) and {105} regions as a function of the nondimensional island 
volume. Prepyramids extending from V = to V = V3 have no {105} regions. T-pyramid solutions occur for V > Vi 
with the amount of {105} increasing monotonically with volume. As the island grows, the transition from prepyramid 
to T-pyramid is energetically favorable for V > V-z, and has no barrier for V > V3. The inverse transition, for a 
shrinking island to change from T-pyramid to prepyramid, is favorable for V < V<i and has no barrier for V <V\. 
Throughout the range V% < V < V3 the transition from prepyramid to T-pyramid is characterized by a jump in the 
{105} width. This behavior is qualitatively similar to the Sid data in Fig. EJb) which shows the abrupt transition in 
the {105} area when the island volume is near 1000 nm 3 . 

If the activation energies are large compared to thermal fluctuations, the transition from prepyramid to T-pyramid 
will occur at V3, while the transition from T-pyramid to prepyramid will occur at V\. Consequently, the prepyramid/T- 
pyramid transition will demonstrate hysteresis behavior as illustrated in Figs-HHb) and (c). Thus, the model suggests 
that the observed difference in the transition size between samples Sid and Shd [see Figs. |5fa) and (b)] can be 
understood if islands in samples Sid are growing, while the small T-pyramids in sample Shd [see Fig. Efd)] are in the 
process of shrinking. This is a probable scenario, since in sample Shd large pyramids act as sink of material and grow 
at the expense of smaller islands. This is analogous to the anomalous ripening process observed by Ross et tdmi in 
the case of the transition from pyramid to dome. Then Figs. Efa) and (b) allow a direct comparison of the volume 
distribution in the case of growing and shrinking T-pyramids in qualitative agreement with the model. The Sid data 
would correspond to the "forward" transition in which T-pyramids occur for V > V3 , while the Shd data sample some 
of the "backward" transition in which T-pyramids can also occur for V\ < V < V^w^ 

Figures EJd) and (e) show the plots of the "area" (width in this two dimensional calculation) as a function of height 
for the theoretical solutions. Here we plot the shape parameters from the theory, taking the transition to occur at V3. 
As with the data, the width of the prepyramid solutions is relatively constant for small island heights. The transition 
to T-pyramids is marked by a discontinuous jump in the height, leading to well-separated curves for prepyramids 
and T-pyramids. Finally, as T-pyramids grow, the area of the {105} regions grow while the area of the (001) region 
shrinks. 

In conclusion, we have investigated by STM the early stages of island formation and faceting during sputter 
deposition of SiGe alloys on Si(001). Besides confirming previous observations obtained with different growth and 
characterization techniques, our study reveals new details of the shape of the pyramid precursors in agreement with 
the predictions of a model recently proposed^ A statistical analysis of the island shape as a function of size allows 
us to follow the morphological transition in which unfacetted prepyramids transform into facetted pyramids and vice 
versa. Our model is able to grasp the main observations reported here. 
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FIG. 1: STM image of a Sio.7Geo.3 layer, 3.8 nm thick, grown on Si(001) at a substrate temperature of 600°C and annealed for 
20 s. Inset: average Fourier transform of six similar images. The observed ring corresponds to an average distance of 125 nm 
between mounds. 




FIG. 2: STM images of Sio.7Geo.3 (a) and Sio.sGeo.s (b) prepyramids and corresponding line scans through the center of the 
islands. Cosine fits of the line scans are also shown. 
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FIG. 3: STM images of Sio.5Geo.s/Si(001) layers 1.5 nm thick, sample Sid (a) and 2.8 nm thick, sample Shd (b), annealed for 
20 s at 600°C. Vertical scale: 8 nm (a) and 10.5 nm (b). The right frames (c) and (d) show the same images, but with gray scale 
according to local surface slope: dark and light areas correspond to {105} facets and orientations close to (001), respectively. 
The arrow in (d) points at a small T- pyramid in the sample Shd- 
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FIG. 4: (a) Magnification of two islands observed in Fig. EJa). Vertical scale: 4.5 nm. (b) Line scan of (a) along the [100] 
direction in correspondance with the black segment on the left side of the image. 
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FIG. 5: Scatter plots for 331 islands (solid squares) on a set of four samples (Sid) obtained by depositing Sio.sGeo.s layers with 
thickness in the range 1.5-1.9 nm on a Si(OOl) substrate at 600° C. The area of the surface occupied by {105} facets (-Ano5)) 
and that of regions with orientation close to (001), (A(ooi)) are plotted vs island volume (a), (b) and height (d), (e). Plots of 
Afooi) and Afros) vs volume for 396 islands contained in sample Shd (see Fig.[3]for details) are also shown in (a) and (b) as open 
circles. Volume and height distributions for samples Su are shown in (c) and (f), respectively. With increasing size, islands 
transform from unfacetted prepyramids (PP) to partially facetted T-pyramids (TP) and then gradually to mature pyramids 
(P). 
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FIG. 6: Schematic representation of the predicted equilibrium shapes (a) and width of (001) and {105} regions vs island volume 
(b), (c) and height (d), (e) from theoretical model. 



